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Abstract 

We consider the class of Markovian processes defined by the equation 
dx/dt = — fix + Zk5{t — tk)- Such processes are encountered in sys- 
tems (hke coalescing systems) where dynamics creates discrete upward 
jumps at random instants and of random height z^. We observe 
that the probability for these processes to remain above their mean 
value during an interval of time T decays as e~^^ defining 9 as the 
persistence exponent. We show that 9 takes the value /3 which thereby 
extends the well known result of the Gaussian noise case to a much 
larger class of non-Gaussian processes. 
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1 Introduction 



Let us consider a stationary random process x{t) and suppose that at t = 
this process is above a given threshold X. What is the probabiUty Q{T) 
that the process has never crossed the threshold X up to time t = T ? In 
many examples of statistical mechanics we observe 

Q{T) ~ e-^^^, T ^ oo (1) 

where 6x is a function of the threshold. The question of determining 9x 
is called the "persistence problem" and has recently been the main interest 
of many papers in the study of out of equilibrium systems. For a short 
review the reader can see Ref. Q and references therein. In many persistence 
problems, only the case where X = (x) is considered and usually one writes 
9 for 6'(a,). 

It often happens that stationary processes can be obtained from physical 
nonstationary systems by rescaling the space variable and taking the new 
time variable t = logr ||l|, where r is the physical time. In terms of the 
latter the persistence probability decays as a power law, thus defining 6x as 
a new independant exponent ^ called the persistence exponent. 

In this paper we present a class of stationary Markov processes x{t) that 
share a universal persistence exponent when X = (x). This class of processes 
is defined by the linear equation of motion 

^ = -Px + rj{t) (2) 

where rj(t) is a noise term consisting of upward jumps uniformly distributed 
along the time axis with density p, 

oo 

77(0= E ^kS{t-tk) (3) 

fc=— oo 

The jump heights Zk are independent identically distributed random vari- 
ables of law R{z) such that {z) exists. 

If the jump heights are all equal to a, the limit a ^ 0, p ^ oo with pa'^ 
fixed transforms the shifted noise r/(t) — pa into Gaussian white noise. In 
this case Eq. (Q) is the usual Langevin equation, which is therefore included 
as a limit in our class of processes. 

The interest of this work is to solve the persistence problem for the pro- 
cess defined by Eq. (^). Simulations have shown that for T large Q{T) decays 
exponentially with an inverse decay rate Ox depending on the threshold and 
on the shape of R{z), such as a Dirac peak, an exponential, or a power law. 
However if X = (x) = p{z) / (3 the result 9 = (3 appears for any law R{z). 
This equality was known, until now, only for the Langevin equation, with 
Gaussian ri{t) [|, |, |. 
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2 A simple method for the calculation of 9 



In this section we show that Ox is the lowest eigenvalue of a differential op- 
erator. Then when the threshold is set to X = (x) , the persistence problem 
is solved straightforwardly. 

2.1 Transformation into an eigenvalue problem 

Let P*{x,t) be the probability distribution of x at time t. It is the solution 
of the master equation 

POO 

dtP{x, t) = (5d^{xP{x, t))+p / P{x - z, t)R{z) dz - pP{x, t) (4) 

Jo 

For the persistence problem only the subclass of functions x{t) that are 
greater than X during the time interval T is relevant. Let P*{x,t)dx be 
the probability that x{t') has not crossed the level X for t' G [0, t] (given 
that x(0) > X) and that x{t) is in the interval [x,a; + dx]. P^{x,t) is 
then solution of the same master equation but with the absorbing boundary 
condition P^{x,t) = when x < X. Because of this condition, P*{x,t) 
will decay to zero as t ^ oo. We will write it as a superposition of decay 
modes P^(x)e~'^^*, where P^{x) is an eigenfunction satisfying the eigenvalue 
equation corresponding to Eq. (^) 

d 

- XxP+ix) = (3—{xP^{x)) +p P^{x- z)R{z) dz - pPlix) 
dx Jo 

with the condition P^{x) = when x < X. Let A^^^ be the eigenvalue of the 
slowest decaying mode. Then we have 

POD 

Q{T)r^e-^xT P^*{x)dx as T ^ oo (6) 
Jx 

so that 

Ox = X*x (7) 

if P^' {x)dx exists. 

Finally let us call P^(x) = P^*(x) the persistence mode. It is necessary 
to have 

POD 

P^ix) > for X > X and / P_^(x)dx < oo (8) 

in order to define P^{x,t) as a probability distribution. Moreover the def- 
inition of the persistence problem ensures that x(0) is chosen within the 
stationary distribution of the process. This distribution is such that (x) 
exists. Therefore we will only focus on solutions -P+(x) with 

POO 

/ xPf(x)dx<oo (9) 
Jx 

As a result solving the persistence problem is now equivalent to calcu- 
lating the lowest eigenvalue of Eq. (|5|) that satisfies Eqs. (^) and (^). 
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2.2 Calculation of 9 

We will now show that A*^.-^ = (3 irrespective of the jump size distribution. 
To this end we define 

fxx= x^PU^)dx, n = 0,l,... (10) 

and we integrate Eq. (|5|) on x from X to oo. As P^{x) = when x < X, we 
obtain the first relation 

Ax/l°i = /3XP^(X) (11) 

We now multiply Eq. (^) by x and integrate as before. After one integration 
by parts we get 

Ax/S = PX^P+{X) + + p)f^^i -p Tdx r dzxP^ix - z)R{z) 

Jx Jo (^2) 

Using the absorbing boundary condition, we can write 

roo roc roo roc 

/ dx dz xPHx - z)R{z) = dyP^{y) dz{y + z)R{ 
Jx Jo Jx Jo 



z] 

(13) 



and the RHS of this equation simply reduces to fxx~'^{^)f\x- Then Eqs. ( p^ ) 



and (ll3|) give the second relation 

Ax/S = PX^PiiX) + - (14) 



We then eliminate P^(X) from Eqs. (p4|) and ( |T^ ) and obtain finally using 
p{z) = I3(x) 

(/S-^/fi)Ax=/3(/S-W/S) (15) 

which is valid for any Ax for which < oo. One can remark that Eqs. (|lT| ) 
and (0), and therefore Eq. (15), could not be obtained if the jump heights 
were not all positive. So this property is essential in the present work. 

Let us now take X equal to the average of the process: X = (x). Eq. ( [T5| ) 
shows that we must have 

A(.>=/3 or /i;V/l(i> = (16) 
If in Eq. (|l6|) we set A^^.) = A*^^ = 9, we see from Eqs and (jl^) that 

ffHx)/ J'eix) average of the process x in the persistence mode, i.e., given 

that it remains always above X = (x). This quantity is therefore larger 
than (x) and the second one of Eqs. ( p!^ cannot be satisfied. Finally the 
only possibility is A*^^ = /3 which proves the main result of this note. 
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For A^^.^ 7^ A*^,^ Eq. (16) shows that 



dx(x-(x))P^(x) = A(,)>A^,) (17) 

{x) 

It follows that for A^^.^ ^ A*^^ any eigenfunction P+(x) with < co must 



have at least one zero for (x) < x < oo. One can also see Eq. (17) as an 
orthogonality relation between two eigenvectors of Eq. ^ that are x — (x) 
(with the eigenvalue —2(3) and P^(x). 

One final remark is necessary. We have studied the persistence prob- 
lem in the most natural situation, that is when the condition /{^],^ < oo is 
fulfilled. However one could be interested in other cases where the initial 
distributions P{x, 0) has a large tail for x — > oo that do not ensure this con- 
dition. If one starts with such a distribution, one can expect the persistence 
problem to be governed by the shape of this initial distribution, with dif- 
ferent answers depending on the choice of the latter. Therefore large initial 
distributions where f'^^ = oo are simply cases of no interest in our study. 

2.3 Perturbation for X close to (x) 

Very few cases in statistical physics offer the possibility of calculating the 
quantity 

(^)x = (18) 

which is the average of the process in the persistence mode. However as 

fe{x) ~ fjsix) assume that 

„(n) _ „(n) fin) /.qn 

Jex — J i3(x) ■lex y^^) 



where f^'^ ^ if |X - (x)| ^ 0. Rewriting Eq. for A^ = A^^ one gets 

f^^^ - (x) f 
Jex ^ Jex 

which can also be written 

= Ljx) (21) 

(x}x-{x) 

Finally one can suppose that X is close to (x) and perturbs Eq. (^) to find 

9x = P{l+ ,^~^_^] , +...) \X-{x)\^0 (22) 

Unfortunately we cannot go further here. We do not know how to calculate 
except in the Gaussian limit which will be the focus of next section. 
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3 The special case of Gaussian noise 



Our goal here is to illustrate the method of Sec. 2 in the special known case 
of Gaussian process. If one considers constant jumps with R{z) = S{z — a) 
and sets pa? = T, one can take the limit a — > with F fixed in Eq. (Q) and 
find 

dtP{x, t) = I3d^{{x - {x))P{x, t)) + ^d,.P{x, t) + 0(a) (23) 

In this limit Eq. ( p3| ) leads to the Fokker-Planck equation for the zero average 
process x{t) — (x) which hence becomes the stationary Gaussian Markovian 
process also called the Ornstein-Uhlenbeck process. This shows that the 
Gaussian case, which is of high interest in statistical physics, is included as 
a special limit in the class of processes studied in the present work. 



3.1 Calculation of 6* 

For simplicity we will now write x instead of x — {x) and X instead of 
X — (x). We thus consider the zero average Gaussian process x{t) and we 
want to determine the persistence exponent 9x related to the probability of 
having never crossed the level X upto time T. As in Sec. ^ it is equivalent 
to find the lowest eigenvalue A^^^ of 

(24) 

The solution of this equation is 

P^{x) = e-y"'\AiD^/p{y) + A2D^/p{-y)) (25) 



where y = \J2(3/Tx is dimensionless and D^j is the parabolic cylinder func- 
tion of index u [^. The integration constants Ai^2 are such that 

AiD^,p{Y) + A2D),ip{-Y) = (26) 



with Y = ^J2I3/TX. After one performs the same integrations as in Sec. 
Eq. ( p^ gives the two relations 



F dP^ 



2 dx 



(27) 

X 



Ax 



(28) 



X 



Eliminating 



X 



from Eqs. (|2^) and (^) one finally gets 



p(0) 



(1) 



(29) 
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which is strictly equivalent to Eq. In the same way as in Sec. 2 it 

now follows that if X = then 9 = \q = P, which is the known result of 
Gaussian persistence |2|, ^. 

Another remark can be made here. The solution of Eq (|2^ ) that vanishes 
at the origin is proportional to (with Ai = —A2 = A) 

P{y) = Ae-y'/\D,/^iy) - D^/p{-y)) (30) 

which can be reduced to 

P{y) = A'ye-y'l^ ^F, ((1 - A//3)/2, 3/2, y^l) (31) 

where A' is simply a numerical constant |^. Using the result Q 

/ e-"n^/2 iFi{{l-u)/2,3/2,u) du = u>l (32) 
Jo 

one has immediately f'^^ = for any A > /?, which is a direct proof of 
Eq. ( p!?!) in the Gaussian case. 



3.2 Perturbation for small X 

Our method offers a faster way than the one presented in 1^] , in the Gaussian 
case, to calculate 9x when X is close to (x) (which is here). 

As P^{x) = A'ye-y'/'^ |], one finds (x)^^.^ = ^ttT/AI3 and Eq. (H) gives 
immediately 

6x = /?(! + + ...), ^ ^ (33) 

which is the well known perturbative result for the Gaussian persistence 
exponent 0]. 



4 Application 

We give here some examples of Markovian but non-Gaussian processes de- 
fined by Eq. (|). 



4.1 ID Percolation 

Let us consider an infinite ID lattice. Each site can be occupied or empty. 
At an initial time, say t = 0, we start with a totally empty lattice. Then, 
at each time step, one site is chosen randomly among non occupied ones 
and is then set occupied forever. The density of occupied sites p is therefore 
an increasing function of time. It goes from to 1 as time goes from to 
00 (in the limit of an infinite lattice, which is the case we consider here). 
We define clusters in a standard way, that is as sets of adjacent occupied 
sites. It is then possible to study percolation on this lattice as a dynamical 
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process, de Freitas and dos Santos Lucena have already performed this 
study numerically. 

The interest here is to consider a semi-infinite lattice and the cluster that 
contains its origin. When p goes from to 1, the size s{p) of this cluster 
will grow from to oo. It is very easy to see that 

{s) = -^ and \/A^=-^ (34) 
1 — p 1 — P 

The evolution equation for s{p) is also quite simple. This random process 
is constant except when new sites are connected to the cluster which creates 
upward jumps of height j. Hence 

ds °° 

T = y2^k^^'P~Pk) (35) 

where k is the jump index. If R{j)Ap is the probability to have a jump j 
when the density of occupied sites goes from p to p + Ap then 

R{j)Ap = pi-Hl-p)-^ (36) 
1-p 

The quantity Ap/{1 — p) is simply the probability to make a jump during 
Ap. The normalised distribution law for the jumps is then 

RU)=I^-\1-P) (37) 

and the density of jumps is p{p) = 1/(1 — p). As p tends to 1, the typical 
time scale of the dynamics of the process s{p), which is ~ 1/p, goes to zero. 
The aim is now to deal with a stationary process. If the new time scale is 

t = log (38) 
1 — p 

the jump times tk are now uniformly distributed with a density 1 along the 
t axis. Then we rescale the process s{p), setting 

.W = -£|, and = ^ (38) 

and we take the limit p — > 1 or equivalently t ^ oo. The jumps z have now 
the probability density R{z) so that R{j)Aj = R{z)Az. In the limit p — > 1 

R{z) = e-' (40) 

If we write Eq. ( ^5| ) in terms of x, t, and z, the evolution equation is finally 

dx °° 

= -xit) + J2 ^k5{t - tk) (41) 

fc=0 

in the limit t ^ oo. This gives a second example of noise of Eq. (^) with an 
exponential distribution law for the jump height. The value of the persis- 
tence exponent is then 9 = 1. 
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4.2 Other examples 

A second example of noise defined by Eq (P) is given by the study of random 
walks. The time t(S) needed for the walker to visit S distinct sites on a 
ID lattice is a random process which, under the same kind of scaling as in 
the last subsection, satisfies Eq. (|2|) ||^]. The jump height distribution R{z) 
is then a power law with a cutoff. In this case the value of ^ is 2. We will 
not go further about this result in the present paper and the reader can see 
references [Q] and P for more details. 

Finally the process x{t) defined by Eq. (H) is simply the electric signal 
that one would observe in an RC circuit receiving random charges due to, 
for instance, the presence of radioactive elements. The persistence exponent 
is then 9 = 1/RC, that is the inverse relaxation time of the circuit. 

5 Conclusion and discussion 

With a very simple calculation we are able to find the persistence exponent 
of a large class of stationary Markov processes which includes the Gaussian 
case. As this quantity is universal within this entire class, it demonstrates 
the existence of universality classes related to non-Gaussian persistence ex- 
ponents. 

To go further one should study the properties of the master equation 
(Q) because the level-dependent exponent 6x is entirely determined by the 
knowledge of {x)-^. However the nontrivial asymptotic behaviour 9x ~ 
XlogX + X^^^{logXf^^ + X ^ oo, found in Ref. [g when the jumps 
are constant shows that this investigation is likely to be very difficult. 
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